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We present a renormalization- group calculation based on model F for the 
' thermal conductivity Xt{T, Q) in the presence of a homogeneous heat current 

, Q and gravity. For temperatures below T\ we obtain a large but finite thermal 

\ conductivity corresponding to superfluid ^He with dissipation. Furthermore, 

we consider the self-organized critical state where the effects of gravity and 
heat current cancel each other so that the distance from criticality AT = 
T(z) — T\{z) is constant in space and a function of Q. We compare our 

■ theoretical results with recent experiments. 

g ! PACS numbers: 67.40.Pm, 64.6O.Ht, 66.60.+a, 64.6O.Lx 
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Liquid ^He at temperatures T near the superfluid transition Tx in the 
^ ' presence of a homogeneous heat current Q offers the possibiHty to study dy- 

■ namical critical phenomena under nonequihbrium conditions. Theoretical^'^ 
and experimental^ investigations in the critical regime have shown that the 
superfluid transition is shifted to lower temperatures by a nonzero heat cur- 
rent Q. The experimentally observed critical temperature Tc{Q) is lower 
than the theoretically predicted T\{Q), i.e. Tc{Q) < T\{Q) < Tx. While for 
T > Tx{Q) the helium is normal fluid and for T < Tc{Q) it is superfluid, Liu 
and Ahlers^ found a new dissipative region for temperatures in the interval 
Tc{Q) < T < Tx{Q). This observation indicates that at a finite heat current 
Q the transition from normal fluid to superfluid helium may happen in two 
steps with two critical temperatures Tx{Q) and Tc{Q). A similar experiment 
was performed by Murphy and Meyer. ^ 

In this letter we present a renormalization-group (RG) theory for the 
thermal conductivity At = Xt(T,Q) and the temperature profile T{z) of 
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liquid ^He in the presence of a nonzero heat current Q and gravity for all 
temperatures above and below T\. We intend to find an explanation of the 
experimental observations. In a previous RG theory^ At(T, Q) has been 
calculated in the normal-fluid region for T > Tx without gravity. Here we 
extend these calculations into the superfiuid region T < Tx. While in the 
previous theories^' ^ the dissipation of the superfluid-normal-fluid counterflow 
related to Q was neglected and hence infinite thermal conductivities were 
obtained for T < Tx{Q), here we present a RG theory for the first time 
which includes the dissipation so that we obtain a large but finite \'t{T,Q) 
fovT<Tx{Q). 

The experiments^"^ have been performed on earth under the influence 
of gravity. While in most cases gravity is negligible, for small heat currents 
below IjLtW/cm^ and temperatures closer than 0.1 juK to Tx gravity effects 
are observable and become important. Via the hydrostatic pressure variation 
induced by gravity the superfiuid transition temperature Tx = Tx{z) depends 
on the altitude coordinate z where the gradient is dTx/dz = +1.273//K/cm.^ 
In the experiments^"^ the heat current Q flows vertically from bottom to top 
to avoid convection. The local thermodynamic state of the helium depends 
on the temperature difference AT(2;) = T{z) — Tx{z) and on Q. Since here 
the gradients VT and VTx have opposite signs, the effects of gravity and 
heat current are additive. 

If the heat current flows downwards from top to bottom then the gradi- 
ents VT and VTx have the same sign. The effects of gravity and heat current 
may cancel each other so that the difference AT = T{z) — Tx{z) is constant 
in space. This configuration has been considered by Onuki^ and suggested 
for an experiment by Ahlers and Liu.^ A self-organized critical state may 
form which is homogeneous over a large region in space. The experiment has 
been realized recently by Moeur et al.^*^ The constant temperature differ- 
ence AT = AT{Q) has been measured as a function of the heat current Q. 
We calculate AT{Q) within our theory and compare with the experiment. 

In liquid ^Hc close to the superfiuid transition heat-transport phenom- 
ena are described by model F^^ which is given by the Langevin equations 
for the order parameter ip{r,t) and the entropy variable m(r,t): 

dil) 5H . ,6H ^ 



where 



H = /'<i''r[i'-o(2)l*P + i|V*p+a„| 
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+5X0 + 7om|V'P - horn] (3) 

is the free energy functional and 9^ and 9jn are Gaussian stochastic forces 
which incorporate the fluctuations. The heat current Q is imposed by bound- 
ary conditions. The gravity is included via the temperature parameter To{z) 
in (3) which is related to Tx(z) and depends linearly on the altitude z. Usu- 
ally, the model is treated by field theoretic means. Critical fluctuations close 
to T\ are taken into account by renormalization and application of the RG 
theory. All the renormalized coupling parameters of model F depending on 
a RG flow parameter and the nonuniversal amplitudes are known^^ so that 
model F can be used for explicit calculations of physical quantities in the 
critical regime without any adjustable parameters. 

The previous field-theoretic approach^' ^ is constructed as a perturba- 
tion theory in terms of Feynman diagrams starting with the mean-field so- 
lutions ipraf{z) and minf(z) of (1) and (2). The quantities ns{z) = (|V'P) and 
3s{z) = (Im [■?/;* V'V']) must be calculated perturbatively so that eventually 
T{z), \t^{z), and hence At(T', Q) can be determined. Below T\{Q) in the 
superfluid region the ansatz Vmf(-2^) = 77 e**^^ and m,^i{z) = const, has been 
used^ which, however, implies a constant temperature proflle. Consequently, 
the gradient VT is zero and the thermal conductivity At is inflnite so that 
dissipation of the supcrfluid-normal-fluid countcrflow is not included. On the 
other hand, in the superfluid state dissipation occurs by creation of vortices 
which implies strong fluctuations of the phase of the complex field V' so that 
(V') = 0. This fact would require the ansatz V'mfC-z^) = which implies that 
mrai{z) is a linear function of z. By using this ansatz the thermal conduc- 
tivity At (7", Q) was calculated previously® for nonzero Q in the normal-fluid 
region where T >T\. Since ns{z) and Js(^;) were evaluated by using the free 
■0-field Green's function, this approach is is invalid below T\. Thus, dissi- 
pation in superfiuid helium is not accessible by perturbation theory starting 
with a mean-field solution of the model-F equations. 

To overcome this problem we suggest the following idea and approxi- 
mation. We generalize model F by replacing the complex field V' by a vector 
* = (V'l) ■ ■ ■ I'tpn) of n complex fields. While model F is recovered for n = 1, 
here we consider the limit n ^ 00. It is well known that models involving 
Ginzburg-Landau functionals like the 0^ model can be solved exactly in this 
limit (see e.g. Ref. 13). The same is true also for model F. Thus, in the 
limit n ^ 00 the complete temperature profile T{z) and the thermal con- 
ductivity At(T, Q) can be calculated exactly. Eventually we set n = 1, apply 
the renormalization-group theory to include critical fluctuations and use the 
result as an approximation.^^ However, first the Green's functions must be 
calculated. While for the m-field Green's function the free Green's function 
is obtained in the limit n — > 00, the V'-field Green's function is determined 
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Fig. 1. The ip-field self energy in the limit n ^ oo. 

by the Dyson equation together with the self energy shown in Fig. 1. In 
the limit n ^ oo only the tadpole diagram survives, where however the full 
propagator line is thick and identified by the exact V'-field Green's function. 
The dashed line represents the interaction mediated by the parameter uq 
in (3) and the m-field Green's function. By solving the Dyson equation we 
obtain self-consistent equations for two effective z-dependent parameters re- 
lated to the temperature profiles T{z) and Tx{z). These equations involve 
ns{z) and Js(2;) as one-loop integrals which are evaluated approximately by 
using the fornuilas of Ref. 6. 

The calculations are very long and complicated and will be published 
elsewhere. 1^ Eventually, we obtain the following results. The temperature 
difference and the thermal conductivity are given by 



AT{z) = T{z) - Tx{z) = TxTa^/\ (1 + ^u[t]A) 
9okB {l + if[r]A-i) 



87rC(r) tj[t]F[t] 



respectively, where 



A 

a 
Ap' 
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-1/6 



C''^-i/2(C) - 1 , 



2^1/ 
1 
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T7?-^i/2(C) - 1 , 

w'\tV^\Ht]w'\t]) \A-/\t]w'\t]) 



.47(7 

87r7[r]F[r] ^mr)?\ 
l + i/[r]^i \gokBTx) ' 



(1 + 8n[r]Ai) , 



T Ta dz 
^1/3 [8 + C-16n[r]AC] = 1 



(4) 
(5) 

(6) 
(7) 
(8) 

(9) 
(10) 

(11) 



Eqs. (6)-(ll) are six equations for seven unknown variables o", r, A, Ai, 
Ap', p'^. Thus, there will be one parameter which can be varied arbitrarily. 
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Preferably we choose for this purpose and determine the other six variables 
by solving the equations. Here u[t], 7[t], ^^''[t], w"[t], F[t], and /[r] are the 
dimensionless renormalized couplings of model F which are known functions 
depending on the renormalization-group flow parameter r.^^ Furthermore 
we need the correlation length ^(r) = CoT""^ with i/ = 0.671 and = 
1.45 X 10~^ cm and the parameter = 2.164 x 10^^ s^^ which is related 
to the entropy density at Tx-^^'^^ The amplitudes A and Ai represent the 
renormalized one- loop contributions by Ug and Jg, respectively. The one- 
loop integrals are represented by the function J-'p{C) which is defined by the 
integral 



This function is known from the previous renormalization-group theory and 
is derived with the assumption that the curvature of the temperature profile 
T{z) can be neglected on the scale of a few correlation lengths ^{t). In a 
renormalization-group theory there is one equation needed which fixes the 
flow parameter r. Here this equation is given by (11). While the flow- 
parameter equation may be arbitrary in general, here (11) is an optimal 
choice. For :§> -|-1 and ( <^ —1 wc recover the standard flow-parameter 
equations of the equilibrium theory^^ for T above and below Tx. 

In (4)-(ll) the space coordinate z does not appear explicitly. Instead 
we have ^ as the variable. After eliminating ( from (4) and (5) we obtain 
the thermal conductivity At = \t{AT,Q) for given heat current Q as a 
function of the temperature difference AT = T{z) — Tx{z). Eventually, the 
temperature profile T{z) is obtained by solving the heat-transport equation 
Q = -AtVT differential equation. 

Now, by solving the equations (4)-(ll) numerically we obtain the fol- 
lowing results. The thermal conductivity At(AT, Q) is shown in Fig. 2 in 
a semi- logarithmic plot versus AT for heat current Q = 100 nW/cm^. The 
solid line is our theoretical prediction while the crosses represent the exper- 
imental data of Day et al.^^ While the equilibrium superfluid transition is 
located at AT = 0, the depressed critical temperature /STx{Q) = Tx{Q) — Tx 
is indicated by the arrow in Fig. 2 and located close to the position where the 
solid line has its maximum slope. Clearly, the agreement between our theory 
and the experiment^^ is quite good for AT > ATx{Q). We note that our the- 
ory does not have any adjustable parameters. However, for AT < ATx{Q) 
the thermal resistivity is very small so that here experimental data are not 
available with sufficient accuracy on a logarithmic scale. Thus, for these low 
temperatures our theoretical curve is a prediction. While the effects of grav- 
ity are included in (10) by the gradient of Tx{z) and the heat current is very 
small, it turns out that the solid line in Fig. 2 is only slightly influenced by 
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Fig. 2. The inverse thermal conductivity Ar^ (AT, Q) (thermal resistivity) 
as a function of ^T{z) = T(z) — Tx(z) for heat current Q = 100 nW/cm^. 
The solid line is the prediction of our theory (4)-(12). The crosses represent 
the experimental data of Day et al.^^ where Q flows upwards. The arrow 
indicates ATx{Q). 

gravity which can be neglected. We find qualitatively similar curves for the 
thermal conductivity At (AT, Q) like this one shown in Fig. 2 also for other 
heat currents Q, where gravity effects are negligible for Q > 100 nW/cm^. 

Liu and Ahlers^ obtained the thermal conductivity At(AT, Q) indi- 
rectly for Q = 42.9 ^W/cm^, which is shown in Fig. 1 of Ref. 4. Again 
for AT > ATx{Q) our theoretical prediction agrees with this result, while 
for AT < ATx{Q) it does not agree. In the dissipative region ATc{Q) < 
AT < ATx{Q) a considerably smaWer thermal conductivity is observed than 
our theory predicts. The discrepancy may possibly be an artifact of the 
experiment caused by a surface effect, because in the experiment"* endplate 
thermometers were used. In a more recent experiment, where only sidewall 
thermometers were used, the dissipative region was not found. A thermal 
conductivity Ax was measured for AT < ATx{Q) which qualitatively has the 
same temperature dependence as our theory but is about 20 times larger. 
Thus, a discrepancy remains, but now the opposite way round. Nevertheless, 
our present theory constitutes a progress and an extension of the previous 
theory to lower temperatures. 
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Fig. 3. The temperature difference AT = T{z) — Tx{z) for the self organized 

critical state where the heat current flows downwards. The sohd line is the 
prediction of our theory (4)-(12). The crosses represent the experimental 
data of Moeur et al.-'^^ 

The self-organized critical state is obtained if the heat current flows 

downwards where Q is negative and the temperature difference AT{z) = AT 
is constant. Thus, the gradient of AT{z) is zero which implies = in (10). 
Solving (4) and (6)-(ll) with the constraint p'l = and eliminating all the 
seven unknown parameters, we eventually obtain AT = AT{Q) as a function 
of Q. The result is shown in Fig. 3. Again, the solid line is our theoreti- 
cal prediction, while the crosses represent the experimental data of Moeur 
et al.^*^ While there are no adjustable parameters, the agreement between 
the theory and the experiment is quantitatively very good for heat currents 
Q < 1.5 //W/cm^. However, for larger heat currents the experimental data 
deviate from our theoretical prediction. For a given AT the larger observed 
Q means that the experimental thermal conductivity At is larger tlian the- 
oretically predicted. 

The strength of the critical fluctuations is indicated by the size of 
the correlation length ,^(r) = £,ot^^ ■ We find that ^(r) is maximum 
for AT = AT\{Q). Furthermore, we have calculated the specific heat 
Cq{AT, Q) at constant heat current. We find a strong and relatively sharp 
maximum at AT = AT\{Q). While the previous theories^'^' predict a sin- 
gularity at T\{Q) for the specific heat, here the anomaly is rounded due to 
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the finite temperature gradient. On the other hand, our theory does not 
predict anything unusual in thermodynamic and transport quantities at the 
lower critical temperature Tc{Q) observed in the experiments.^' ^'^'^ 

Furthermore, we have calculated the temperature profile T{z) of the 
superfluid-normal-fluid interface. For heat currents Q > 1 fxW/cm^ we ex- 
pect that the effects of gravity are small. In the normal-fluid region wc find 
that for these heat currents gravity is indeed negligible. However, in the su- 
perfiuid region gravity effects are important even for larger heat currents. In 
microgravity we predict that T{z) decreases monotonically without a lower 
bound for increasing z which is due to the fact that the thermal conductivity 
At remains finite. In gravity where the heat current flows upwards we find 
a temperature T{z) = Tc{Q) where the gradients VT{z) and VTx{z) have 
equal magnitude but opposite signs. For Tc{Q) < T{z) < Tx{Q) the gra- 
dient VT{z) dominates and in this region the temperature T{z) decreases 
monotonically as in microgravity. However, for T[z) < Tc{Q) the gradient 
VTx{z) dominates and for increasing z the temperature approaches the lim- 
iting value lim2:_»oo ^(2;) = Too{Q) which is not much below Tc{Q). Thus, 
our theory indeed yields a dissipativc region and a superfluid region nearly 
without dissipation separated by a temperature Tc{Q). The obtained shift 
ATc{Q) is nearly the same as AT{Q) in the self-organized critical state, 
which has been found also in the experiments.^' ^'^^ Our theory predicts 
that Tc{Q) is not a property of liquid ^He but it is an artifact which is gen- 
erated by gravity and occurs when integrating the heat-transport equation 
to obtain the temperature profile. In microgravity Tc{Q) will not be present. 

From these observations we conclude that our theory concerning 
Tc{Q) agrees qualitatively with the experimental findings.'^"''''^' However, 
there are two major quantitative disagreements. First of all we predict 
ATc{Q) ~ —Q^'^^ for Q > 1 //W/cm^ which disagrees with the experi- 
mentally observed^ Arc^exp(Q) ~ —Q^'^^. Secondly, in the dissipative region 
Tc{Q) < T{z) < T\{Q) the thermal conductivity At is much larger than 
found in the experiment^ which implies that the spatial extent Az of the 
dissipative region is much larger. We find Az > 0.25 cm for Q > 10 jiW /cm' 
which is about the sample size in the experiments.^"^ However, it is unclear 
if the dissipative regions and the related definitions of Tc{Q) in theory and 
in the experiments have the same origin. The region of enhanced dissipa- 
tion may possibly be a surface effect and not a property of the bulk helium, 
because it was not observed in the recent experiments^ using only sidewall 
thermometers. Thus, it is unclear if the discrepancies are due to the theory 
or due to the experiment. 

For the understanding of the nature of Tc{Q) and of the dissipative re- 
gion further experimental and theoretical investigations are necessary. The 
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following questions should be addressed in future experiments. Does Tc{Q) 
depend on gravity? Is Tc{Q) absent in microgravity as our theory predicts? 
How large is the spatial extent Az of the dissipative region? (For this pur- 
pose, the spatial dependence of T{z) should be investigated carefully.) Do 
thermodynamic quantities show any anomalies at Tc{Q) or not? On the 
other hand, the discrepancies may be caused by the crucial assumption of 
our theory: we have generalized model F by introducing an n-component 
complex order parameter and obtained our approximation by taking the 
limit n — > CO. It turns out that for n > 2 already in mean-field approxi- 
mation the supcrfluid-normal-fluid counterflow is unstable, while for n = 1 
(model F) it is metastablc. This fact may cause a qualitative difference 
between model F and our approximation. 

Nevertheless, our approach appears to describe the essential part of the 
physics correctly: we find dissipation for the superfluid-normal-fluid counter- 
flow and obtain a finite thermal conductivity At in the superfluid region for 
nonzero heat currents Q. Since dissipation in the superfluid state is necessar- 
ily connected with the creation of vortices, our approximation must include 
vortices in an indirect way. We have applied our approach also to rotating 
helium.^^ In this case, the number of vortices is precisely known and related 
to the rotation frequency, so that the dissipative effect of a single vortex line 
can be extracted which is related to the Vinen coefficient^^ B. We obtain 
the simple result^^ B = (4m4//i)r'[r] where T'[t] is a rcnormalized coupling 
of model F depending on the RG flow parameter r. This result has a simi- 
lar temperature dependence and the same order of magnitude as the Vinen 
coefficient B calculated within the renormalized mean-field theory by con- 
sidering a single vortex line.^^ Thus, the large discrepancies of the thermal 
conductivities by about a factor of 20, which we have discussed above, are 
not caused by an incorrect treatment of the dissipative effect of the vortex 
lines. Rather, the discrepancies must be due to the density of the vortices in 
the helium, which in the homogeneous heat flow appears to be a very scnsity 
quantity possibly influenced strongly by the kind of the approximation and 
by the conditions of the experiment. 

Finally, our approach can be tested in thermal equilibrium for Q = 
by calculating e.g. the speciflc heat Cp. Our result for Cp agrees with the 
results of the previous theories^^ within the related accuracies. While the 
correct critical behavior is incorporated by the RG theory, we obtain the 
amplitude function F[u] of the speciflc heat, which below Tx agrees with 
the amplitude function F^[u] of Ref. 12 in renormalized mean-fleld theory 
(zero-loop approximation) and above Tx agrees with -F+fu] of Ref. 12 in one- 
loop approximation. Thus, we conclude that our approach is valid at least 
in thermal equilibrium. 
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